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Abstract 

Radio-loud quasars usually contain parsec-scale nonthermal jets. The most compact 
emission region (“the core”), and perhaps some of the moving “knots,” are expected to 
be efficient producers of inverse Compton scattered X-rays and 7 -rays since many of the 
synchrotron photons will upscatter before escaping. Through multifrequency flux density 
observations and VLBI measurements of angular sizes, one can predict the flux density of this 
self-Compton high-energy emission. It is not always the case that the brightest synchrotron 
sources are also the brightest X-ray and 7 -ray sources. Perhaps a better predictor of 
high-energy brightness is the ratio of hard X-ray to high-frequency laclio emission. 

Using the synchrotron self-Compton relations, we predict the 7 -ray fluxes of several 
sources we expect to be detected by EGRET. More accurate predictions will be made when 
we complete a program of contemporaneous radio-submillimeter and X-ray observations 
during the course of the EGRET all-sky survey. 


1. INTRODUCTION 


VLBI observations of the parsec-scale structure of quasars show that the radio jets 
are extremely luminous, suggesting high energy densities of soft photons. It is expected 
that the synchrotron photons will scatter ofT the relativistic electrons that produced them to 
create X-rays and possibly soft 7 -rays. Second order scattering should produce hard 7 -rays. 
The synchrotron self-Compton (SSC) process has been discussed in the context of quasars 
by several authors (e.g., Jones, O’Dell and Stein , 1974; Ghisselini, Maraschi, and Treves, 
1985; Band and Grindlay, 1986; Marscher 1987). The large X-ray data base provided by 
the Einstein Observatory has made it feasible to test the relationship between radio and 
high-energy emission, at least in a statistical sense. Several authois show that a stiong 
correlation exists between X-ray and radio emission (Owen, Helfand, and Spangler 1981; 
Worrall et al. 1987; Kembhavi, Feigelson, and Singh 1986; Browne and Murphy 1987). 
Also, 3C273 and 3C279, both strong radio sources, have been detected in 7 -rays . These 
results would suggest a one-to-one correspondence between the high-energy flux and radio 
flux. However, a more detailed look at self-Compton theory shows that the X-ray flux density 
F vx will have, in addition to a dependence on the radio flux density, a strong dependence 
on several other parameters (e.g., angular size). Similar calculations can be extended to the 
second order Compton flux densities. The dependence on several parameters is even more 
extreme in this case. Modest uncertainties in these parameters are typical; thus, predictions 
for the 7 -ray flux density made in this manner are quite rough. There is then sufficient 
reason to seek an alternative method of predicting the 7 -ray flux density. We show that if 
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the X-ray flux density is measured, then the ratio F vx IF vr can be used to predict F V 1 /F VX . 
This method was discussed by Jones (1979); we present more precise formulae. Flux densities 
at energies within the EGRET band can then be predicted. 


2. FIRST AND SECOND ORDER SCATTERING RELATIONSHIPS 


For homogeneous sources with simple geometries (e.g., a sphere), the predicted X-ray 
flux density can be calculated using the standard expressions for first order scattering 
emissivities in the Thomson limit (Blumenthal and Gould 1970; Gould 1979; Rybicki and 
Lightman 1979; Hughes and Miller 1991). The first order flux density then has the following 
dependence on the physical parameters of the source: 

F v xc oc N 0 RF„ s ln{ v u /v m ) oc R A D~ 2 ln(u u /v m ). (1) 

Here, R is the radius of the source, and D is the distance. Using expressions for 
the synchrotron flux density and optical depth at spectral turnover for a self-absorbed 
synchrotron source, one can solve for B, the magnetic field stength, and the parameter Ao, 
where the number density of relativistic electrons of a given energy follows the relationship 
N(E) oc NoE~ p , in terms of observable radio source properties. These expressions for the 
physical parameters are then substituted back into the derived equation lor X-ray flux density 
(see Marscher 1987 for details). The major result is that the first order flux density at photon 
energy E^tV keV is highly dependent on these parameters: 

F w x ( 2 ) 


Here, F m and u m are the flux density and frequency at spectral turnover, z is the redshift, 6 
is the Doppler beaming factor (a correction for bulk relativistic motion within the source), 
and the optically thin log/q, 5 vs. log v spectrum is assumed to be a power law of slope 
(1 —p)j 2 between frequencies v m and v u . In the remaining proportionalities, we will ignore 
the slowly varying logarithmic term. 

The second order flux density is then even more strongly dependent on the observable 
parameters (due to an extra factor of NqR): 
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Usually, real sources cannot be well modeled as a single homogeneous sphere; however, 
the VLBI components collectively can be approximated as a composite of homogeneous 
spheres (Marscher 1987). The predicted flux densities for each component are then added 
together and compared to the measured flux density. Such calculations are usually assumed 
to be in the Thomson limit (h»/j n j|< 0 / << me 2 / 7 ) or in the extreme Klein-Nishina limit 
(hvinitiai > mc 2 /~f) because the full expression for the electron differential cross section that 
goes into calculating the emissivities is greatly simplified in these two limits. However, to 



a first approximation, the basic proportionalities of the first and second order flux densities 
on these radio source parameters are independent of these assumptions. More precise 
expressions, which explicitly use the full Klein-Nishina cross section, are discussed below. 

Ideally, we would like to obtain an expression for predicting the 7 -ray flux density that 
is not so highly dependent on parameters that are poorly known. We can see from expressions 
( 1 ) and (3) that the ratio of second order ( 7 -ray) to first order (X-ray) self-Compton 
flux densities equals the ratio of the first order self-Compton to synchrotron flux density. 
Therefore, if the ratio F vx jF vr is measured and assumed to be determined by the SSC 
process, F V1 IF VX can be predicted without knowledge of the radio source parameters. For 
the first order Compton emissivity, it is safe to assume that the Thomson limit applies. 
However, for the second order emissivity the Thomson limit will only apply over a fraction 
of the spectrum. Thus we need to integrate numerically the full expression of Blumenthal 
and Gould (their eq. 2.75). Performing this calculation, we obtain: 
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Here, v\ is the first order photon frequency, is the second order photon frequency, and 
F\{p) is that of Hughes and Miller (1991). The parameter S represents the integral: 

v / dv\ /„(!/, )v\ 2 Jo dq /(</, J'l, J' 2 ) 

M* / lC'2,<7) = IFTT ' y 


Here, is the first order self-Compton intensity, K\c = h{>' i)"i 2 is tlie f ,rst order 

intensity without the frequency dependence, / is the complicated function in equation (2.75) 
of Blumenthal and Gould, and q = « T , l(1 J£, l/mc gj • The integral ,9 is only weakly dependent 
on the upper and lower limits of the first order Compton spectrum and spectral index, and 
is of order unity. We can now solve for Fi/^j /Vx hi terms of P vr j F l/r ■ 
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We must emphasize that using this equation does not eliminate all of the uncertainty of 
expression (3). For instance, a rigorous analysis would require calculating this ratio for each 
separate component of the radio source, which is not possible due to the comparatively low 
resolution of current X-ray and 7 -ray instrumentation. Possible source variability in either 
the radio or X-rays will lead to inaccurate estimates of P], r /Pvr unless the measurements 
are nearly simultaneous. In addition, the Thomson limit is not valid over much of the 7 -ray 
range. In this case, numerical calculations are required (see Marscher and Bloom, these 

proceedings). 

We present the results of a numerical computation of the self-Compton spectrum of a 
model uniform source in Figure 1. Note that there is substantial curvature 1 o the spectrum 
of the self-Compton emission, contrary to the popular myth t hat it is a power law over most 

X-ray and 7 -ray energies. 
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Figure 1 vF v spectrum of the synchrotron (left-most curve), first order self-Compton (middle curve), and 
second order self-Compton (right-most curve) emission from a model spherical, uniform, nonthermal source. 
The lower panel gives the slope of the vF v synchrotron (left-hand) and first + second order self-Compton 
(right-hand) spectrum as a function of frequency. 

Equation (6) is still practical, especially if the detailed observations required for analysis 
of expression (3) are uncertain or non-existent. The best procedure would be to use both 
methods whenever the observations allow. 


342 





3. PREDICTIONS 


Using equation ( 6 ) and taking the accompanying caveats to heart, we can predict which 
sources are likely to be detected by EGRET. We would like to find a sample of sources with 
recent simultaneous radio core and X-ray flux density measurements. Such a sample does not 
yet exist; however, Bloom and Marscher (1991) have compiled the X-ray data for all quasars 
and active galactic nuclei observed with both Einstein Observatory and with VLBI We use 
these data as a rough guideline to predict F vy , although these calcuations are only valid 
for the epoch of the X-ray and radio observations, 1979-198E Table 1 shows sources with 
highest values of predicted 7 -ray flux that are above 5 x 10 photons cm " s , the flux 
limit of EGRET (Kanbach 1989). The second column gives the X-ray flux density at 2 keV 
as measured by Einstein ; the third column gives the measured ratio of X-ray flux density at 2 
keV to the core radio flux density at 10.7 GHz in 1980; the fourth column gives the predicted 
ratio of 7 -ray flux density at 100 MeV to X-ray flux density; and the last column gives the 
predicted photon flux in the EGRET band, integrated numerically between 100 Me\ and JU 
GeV . We have included effects due to the Klein-Nishina electron cross-section, which 

is significant at this energy range. Wc notice that some of the sources with high predicted 
7 -ray fluxes have rather modest X-ray and radio flux densities, implying again that it us not 
sufficient to use radio or X-ray flux alone as an indicator of which sources will be brightest in 
7 -rays. We see from equation (4) that the main requirement is that the source be compact. 
Alternatively, from expressions (1) and (3) we see that the 7~ ia y fb' x ' s piopoitiona. to 1 0 
and is therefore strongly dependent on the relativistic electron density within the source. 


Table 1. Predicted 7 -ray Fluxes, Epoch 1980 


Source 

F VI (ft- ly) 

(10- 6 ) 

r v r 

TT 

1 

O 

4 $ 

F~, (10 7 [fliot cm 2 s 1 ) 

NRAO 140 

1.65 

2.8 

3.6 

0.9 

3C 111 

1.87 

1.7 

2.2 

0.6 

3C 120 

5.2 

5.8 

7.7 

6.0 

3C 249.1 

0.47 

13.4 

17.7 

1.3 

ON 231 

0.31 

0.4 

0.5 

0.02 

3C 273 

7.9 

0.8 

1.04 

1.3 

3C 279 

0.40 

0.1 

0.1 

0.01 

3C 334 

0.15 

2.8 

3.7 

0.08 

Mkn 501 

7.4 

11.2 

14.8 

17 

1721+343 

1.4 

14.0 

18.5 

3.1 

3C 390.3 

1.0 

2.4 

3.2 

0.5 

BL Lac 

1.87 

1.9 

2.5 

0.7 

3C 446 

0.48 

0.2 

0.3 

0.02 
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We note that the predicted 7-ray fluxes of several of the sources in Table 1 (especially 
Mkn 501) were above the detection level of COS-B. This result would imply that some 
of these sources (such as the BL Lac objects) are probably creating X-rays through some 
process other than the SSC mechanism. The predicted flux for 3C273, on the other hand, 
is approximately the same as that measured by COS-B (Swanenburg ft al. 1978; see also 
Jones 1979). The predicted flux for 3C 279 is well below that recently measured by EGRET. 
This low prediction, which is based on X-ray and radio flux densities from 1980, is consistent 
with the non-detection of 3C 279 by COS-B. In addition, the radio flux at centimeter and 
millimeter wavelengths was considerably lower at the time of COS-B observations. For 
instance, the total flux density at 90 GHz went from 7-8 Jy in 1976-1978 (Landau et al. 
1980) to approximately 14 Jy in July 1991 (extrapolated from the 150 GHz flux density; 
Robson 1991). 

During 1991, 3C 279 has been flaring at millimeter and infrared wavelengths (Robson 
1991). The self-Compton process acts as an amplifier of the synchrotron flare, leading to a 
very high 7-ray flux, as observed. Other quasars currently undergoing synchrotron flares are 
the bright radio quasar 1C 39.25 and the high F VI jF vr quasar NRAO 110. We predict that 
these sources will also be detected by EGRET. 

4. CONCLUSIONS 

Compton theory for homogeneous sources shows that the second order flux density can 
be predicted if the ratio F ul jF vr is measured. Equation (6) can then be used to estimate 
the 7-ray flux. However, it is clear Irom the predictions in section 1 that a more accurate 
analysis, looking at each component of the radio source and applying the first and second 
order scattering relationships, is necessary to test the dominance of the SSC process. Future 
work that will include multi wavelength VLBI radio and submillimeter radio observations at 
times close to ROSA! observations will allow for more reliable measurements of F vr jF vr . 
and a detailed component-by-component analysis for a number of quasars. Accurate fluxes 
can then be predicted and compared with EGRET observations. 
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